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ABSTRACT 


Analysis  of  data  obtained  during  the  High  Altitude  Clear  Air  Turbulence 
Program  in  the  altitude  range  45 >000  to  70,000  ft.  indicated  that  probability 
densities  and  exceedances  of  root -mean-square  (RMS)  gust  velocities  varied 
significantly  with  the  character  of  underlying  topography  and  to  a  lesser 
degree  with  altitude.  The  relative  frequency  of  occurrence  of  large  RMS 
values  increased  with  terrain  roughness.  The  components  of  the  RMS  gust 
velocities  were  linearly  related, with  the  lateral  and  longitudinal  component 
values  being  larger  than  the  vertical. 

Lengths  of  the  turbulent  regions  ranged  from  less  than  one  mile  to  130  miles, 
with  the  frequency  of  occurrence  decreasing  rapidly  with  increasing  length. 
Because  the  long  regions  were  small  in  number,  relatively  few  power  spectral 
density  curves  were  derived  over  the  frequency  range  necessary  in  computing 
scale  lengths.  Large  scale  lengths  (generally  8,000  to  l6,000  ft)  and 
relatively  steep  slopes  (-1.5  to  -1.7)  were  suggested  for  mountain-related 
cases  and  for  the  lateral  and  longitudinal  components  from  the  few  Bpectra 
available.  Scale  lengths  were  equally  as  large  but  slopes  significantly 
less  for  the  flat  terrain  spectra  and  for  the  vertical  component.  No 
obvious  relationship  existed  between  spectral  shape  and  altitude  or  wind 
direction. 

The  frequency  of  turbulence  varied  from  2.7$  to  4.6$  for  flights  over  flat 
terrain  and  high  mountains,  respectively.  Above  60, 000  ft.  moderate  turbu¬ 
lence  was  not  uncommon  over  mountains  but  rarely  was  noted  over  low  relief 
or  water. 

The  fairly  high  number  of  turbulence  cases  above  thunderstorm  tops  was 
largely  influenced  by  location  of  base  sites  and  by  turbulence  searching 
techniques.  Consequently,  thunderstorm  turbulence  was  not  classified 
separately. 

These  analyses  resulted  in  suggested  changes  in  the  standard  deviations  of 
the  root -mean-square  gust  velocities,  the  ratios  of  turbulent  flight  miles 
to  total  flight  miles  and  the  form  of  the  power  spectral  density  curves  in 
current  use  for  the  design  and  operation  of  aircraft  in  the  45,000  to  65,000 
ft  altitude  range. 
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DPTHCDUCTIOB 

This  report  consists  of  a  discussion  of  those  characteristics  of  high 
altitude  clear  air  turbulence  that  are  of  direct  Interest  to  those  engaged 
In  the  design  and  operation  of  advanced  aircraft.  Specifically,  the  prob¬ 
ability  density  distribution  of  the  root -mean -square  gust  velocities,  the 
proportion  of  flight  miles  In  turbulence,  and  the  scale  of  turbulence  are 
given  as  functions  of  altitude  and  underlying  topography. 

Factors  such  as  the  "scale  length  of  the  turbulence",  distribution  of  gust 
velocities,  turbulence  above  thunderstorm  clouds  and  the  relation  between 
clear  air  turbulence  and  the  flight  direction  with  respect  to  the  wind  and 
mountain  ridges  are  discussed  In  detail. 

The  results  presented  In  this  report  were  obtained  through  the  use  of 
statistical  and  theoretical  analyses  of  the  complete  set  of  data  from  the 
HTCAT  flight  program.  In  many  Instances  the  data  published  by  Crooks,  et 
al.  (1,2)  were  supplemented  with  unpublished  data* 

The  general  organization  of  this  report  Is  as  follovB.  First,  the  prob¬ 
ability  density  distributions  at  the  root-mean -square  gust  velocities 
obtained  from  truncated  power  spectral  density  curves  were  determined  for 
the  entire  HICAT  sample  and  for  the  saaiple  divided  Into  altitude  and  topo¬ 
graphic  categories.  Sext,  ■'.he  shape,  slope,  and  associated  turbulence 
"scale  lengths”  of  the  power  spectral  density  curves  were  determined. 

Once  the  appropriate  mathematical  representations  of  the  spectra  were 
obtained,  the  relationship  was  established  between  the  total  root-mean- 
square  gust  velocities  and  the  root -mean-square  gust  velocities  obtained 
from  the  truncated  spectra.  This  led  directly  to  the  probability  density 
distributions  of  the  total  root -mean-square  gust  velocities.  The  standard 
deviations  of  these  velocities  gives  the  "b’s"  for  the  equations 
caoaonly  used  by  structural  engineers. .  Then  the  proportion  of  flight  time 
In  turbulence  (the  "P's"  of  the  equations  used  by  the  engineers)  are  given. 
The  remainder  of  the  report  discusses  other  characteristics  of  high  alitude 
turbulence. 
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SECTION  II 


PROBABILITY  DENSITIES  OF  ROOT-MEAN -SQUARE  GUST  VELOCITIES 
Types  of  Root -Mean-Square  Gust  Velocities 

Two  different  values  of  the  root -mean-square  gust  velocities  have  been  found 
to  be  useful  for  the  design  and  operation  of  aircraft.  One  of  these  values 
corresponds  to  the  square  root  of  the  area  under  the  power  spectral  density 
curve  between  wavelength  (frequency)  limits  that  are  within  the  range  of  the 
measurements.  The  second  value  of  the  root -mean -square  gust  velocity  fre¬ 
quently  used  is  given  by  the  square  root  of  the  area  under  the  complete 
spectrum.  This  is  the  value  used  in  most  analytical  representations  of  the 
gust  velocities.  Procedures  for  estimating  the  total  areas  if  the  areas 
under  the  truncated  curves  are  known  sire  discussed  in  Section  IV. 

The  HICAT  Sample 

In  the  HICAT  program  the  turbulent  regions  varied  in  length  from  approximately 
1  to  115  nautical  miles.  The  distribution  of  the  lengths  of  the  turbulent 
regions  is  discussed  in  detail  by  Ashbum  et  al.  (3)  and  Ashbum  (4).  The 
procedure  that  was  adopted  for  computing  the  spectra  required  that  the  largest 
wavelength  in  the  spectra  be  represented  by  at  least  twenty  samples.  Thus, 
the  long  wave  limits  of  the  computed  spectra  were  determined  by  the  observed 
lengths  of  the  turbulent  regions.  The  lower  wavelength  limit  of  appro tiraately 
130  ft.  was  determined  largely  by  the  response  of  the  instruments,  the  speed 
of  the  aircraft  and  needs  of  the  designers  of  aircraft.  The  root-mean-square 
gust  velocities  were  computed  from  the  spectra  with  the  long  wavelength  limits 
set  at  1,000  ft.,  2,000  ft.,  4,000  ft.,  10,000  ft.,  and  20,000  ft.  for  all 
turbulent  regions  whose  length  was  at  least  20  times  the  long  wavelength  limit. 

The  size  of  the  cample  of  root-mean-square  gust  velocities  decreased  signifi¬ 
cantly  as  the  upper  truncation  line  was  moved  to  longer  wavelengths.  This  is 
shown  in  Table  I. 


TABLE  I 

Number  of  RMS  Gust  Velocities  Computed  for  Various  Truncation  Limits 

1 


Upper  Wavelength 

Truncation  Limit 

Number  of  Cases 

2,000  ft.  [  RMS  (X  =  2,000  ft)  ] 

197 

4,000  ft.  [  RMS  (X  =  4,000  ft)  ] 

116 

10,000  ft.  [  RMS  (X  =10,000  ft)  ] 

46 

20,000  ft.  [  RMS  (X  =20,000  ft)  ] 

9 

2 


The  criteria  for  each  of  the  cases  listed  in 


le  I  were: 


a)  The  RMS  (X  =  2,000  ft)  gust  velocity  for  the  vertical  component  was 
equal  to  or  greater  than  0.5  ft/sec. 

b)  Reliable-  spectra  for  all  three  gust  velocity  components  were  avail¬ 
able. 

The  division  by  season,  type  of  underlying  terrain  and  altitude  of  the  197 
cases  of  RMS  (X  =  2,000  ft)  is  given  in  Table  II. 


TABLE  II 

Distribution  of  RMS  (X  =  2,000  ft)  Cases  by  Season,  Terrain  and  Altitude 


Winter 

Spring 

Summer 

Autumn 

All  Seasons 

Water 

2 

11 

18 

0 

31 

Flatland 

19 

25 

5 

4 

53 

Low  Mountains 

17 

4 

0 

13 

34 

High  Mountains 

37 

14 

0 

28 

79 

All  Terrain 

75 

54 

23 

^5 

197 

45,000  to  49,900  ft 

6 

5 

4 

6 

21 

50,000  «o  54,900  ft 

31*- 

15 

15 

6 

70 

55,000  to  59,900  ft 

16 

?6 

3 

17 

62 

>  60,000  ft 

17 

8 

0 

16 

41 

Selection  of  Sample 

The  distribution  of  lengths  of  the  turbulent  regions  encountered  in  the  HICAT 
flight  program  (Ashburn,  4)  was  such  that  approximately  85$  of  the  turbulent 
regions  were  long  enough  to  permit  the  computation  of  RMS  (X  =  2,000  ft)  under 
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the  criteria  established  by  Crooks  et  al.,(l,  2).  The  corresponding  percent¬ 
ages  for  RM3  (\  3  10,000  ft)  and  RMS  (X  =  20,000  ft),  where  there  b*d  been  no 
malfunctions  of  instruments,  were  approximately  20^  and  5$,  respectively. 

These  data,  along  with  the  results  presented  in  Tables  I  and  II,  indicate  that 
if  a  statistical  summary  of  the  RMS  gust  velocities  is  to  be  made  by  categories 
of  underlying  terrain,  season  or  altitude,  then  only  the  RMS  (X  3  2,000  ft) 
values  provide  an  adequate  sample  size.  Selecting  the  RMS  (X  *  2,000  ft)  values 
for  the  analysis  has  some  obvious  disadvantages  from  the  point  of  view  of  study¬ 
ing  important  characteristics  of  the  spectra  but  this  choice  had  to  be  made  if 
the  root-mean-square  gust  velocities  were  to  be  investigated  as  functions  of 
altitude,  season  and  terrain. 

Relationship  Between  Lateral,  Longitudinal  and  Vertical  HMB  (X  -  2,000  ft)  Gust 
Velocities 

Correlation  coefficients,  r,  between  the  lateral -vertical,  longitudinal -vertical, 
and  the  lateral -longitudinal  RMS  (X  3  2,000  ft)  values  were  computed  using  the 
following  relationship: 

n  E  x  y  -  J  x  E  y 
i  yi  i  yi 

r  =  ~  . . -  . . — . . . 

f/[n  E  Xj2  -  (  £  x1)2]  [n  2  y^  -  (£  y^2] 


where  the  emanations  are  taken  from  i  =  1  to  i  =  n  and  the  x,  and  y  are  the 
appropriate  RMS  (X  3  2,000  ft)  values.  The  results  given  inthe  second  column 
of  Table  III  were  obtained  by  giving  equal  weight  to  the  values  published  by 
Crooks  et  al  (l,  2).  The  results  given  in  the  first  column  were  obtained  by 
weighting  each  value  according  to  the  length  of  the  time  in  turbulence.  For 
example,  a  turbulent  region  that  lasted  for  500  seconds  of  flight  time  was 
given  five  times  the  weight  of  a  turbulent  region  that  was  traversed  in  100 
seconds.  If  the  turbulence  that  lasted  500  seconds  was  stationary  in  char¬ 
acter,  then  it  could  have  been  divided  into  five  100  second  turbulent  regions. 
This  is  the  Justification  for  weighting  the  RMS  (X  3  2,000  ft)  values  by  the 
time  in  turbulence.  In  general  there  was  no  statistically  significant 
difference  between  the  results  obtained  by  no  weighting  and  weighting.  The 
relatively  high  correlation  coefficients  given  in  Table  III  suggest  that  the 
components  of  RMS  (X  3  2,000  ft)  are  linearly  related  and  hence  may  be 
represented  by 


RMS  (X  3  2,000  ft)L  3  m^  RMS  (X  3  2,000  ft)y  +  81 

RMS  (X  3  2,000  ft)p  =  mg  RMS  (X  3  2,000  ft)y  +  Bg  (l) 

RMS  (X  3  2,000  ft)L  =  RMS  (X  3  2,000  ft )F 
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;ht  Altitude 


TABLE  III 


Correlatioo  Coefficients  Between  Components  of  Gust  Velocity 


Lateral -Vertical 
Water  Longitudinal -Vertical 
.  Lateral -Longitudinal 


Lateral -Vertical 

Flatland  Longitudinal -Vertical 
Lateral -Longitudinal 


Low  Lateral -Vertical 
Mountains  Longitudinal -Vertical 
Lateral -Longitudinal 


High  Lateral -Vertical 
Mountains  Longitudinal -Vertical 
Lateral -Longitudinal 


Lateral -Vertical 

All  Cases  Longitudinal -Vertical 
Lateral -Longitudinal 


45,000  to  Lateral -Vertical 

49, 900  Ft  Longitudinal -Vertical 

Lateral  -Longitudinal 


50,000  to  Lateral -Vertical 
54, 900  Ft  Longitudinal -Vertical 
Lateral -Longitudinal 


55,000  to  Lateral -Vertical 

59,  900  Ft  Longitudinal -Vertical 

Lateral -Longitudinal 


60,000  Ft  Lateral -Vertical 

and  Above  Longitudinal -Vertical 
lateral -Longitudinal 


Correlation  Coefficients 

Weighted  by  Time 

In  Turbulence 

Ron -Weighted 

.899 

.938 

•796 

•835 

•793 

.847 

•965 

-937 

•  942 

.908 

•  946 

-SDL9 

.870 

.845 

.883 

.890 

.896 

.870 

.966 

.960 

.962 

.962 

.964 

.960 

•955 

.950 

•953 

.950 

.950 

.944 

.938 

•  950 

•  948 

.961 

.867 

.881 

.963 

.963 

.964 

.963 

.947 

.947 

•967 

•  943 

.958 

.946 

.958 

.940 

•951 

.951 

•955 

.954 

•971 

•  971 

The  values  of  m^,  m^,  m_,  8_ ,  9^  and  8_  are  given  in  Table  IV.  These  values 
nay  be  used  for  two  purposes.  First,  they  relate  to  the  problem  of  estimating 
the  degree  of  isotropy  of  the  turbulence.  Secondly,  if  the  existence  of  tur¬ 
bulence  is  defined  in  terms  of  the  vertical  component  of  the  EMS  (\  =  2,000  ft) 
then  equations  (l)  may  be  used  to  establish  the  definition  of  turbulence  in 
terms  of  the  lateral  and  longitudinal  components. 

The  data  presented  in  Table  IV  supplement  the  detailed  discussion  of  the  ratios 
of  the  components  of  the  root -mean-square  gust  velocities  that  was  given  by 
Crooks  et  al.  (l)  by  presenting  an  analysis  of  a  larger  sample  and  by  dividing 
the  sample  into  altitude  and  terrain  categories.  This  new  analysis  indicates 
that  the  ratio  of  any  two  of  the  components  of  RMS  (X  =  2,000  ft)  is  not  a 
constant  but  varies  with  the  magnitude  of  the  root-mean-square  gust  velocities. 
Further,  the  analysis  indicates  that  the  relationship  between  any  two  of  the 
components  changes  significantly  with  altitude  and  terrain  category. 

The  values  given  in  Table  IV  may  be  looked  upon  as  useful  in  providing  the 
best  estimate  of  the  relative  magnitudes  of  the  components  of  RMS  (X  =  2,000 
ft).  Tlie  relationships  between  the  components  of  &  are  given  in  Section  III. 

Probability  Densities  of  the  Root -Mean -Square  Gust  Velocities 

The  probability  densities  of  the  root -mean -square  gust  velocities  obtained 
from  the  HICAT  program  may  be  represented  by  histograms  or  by  analytical 
functions.  The  histograms  represent  the  data  more  directly.  Analytical 
functions  are  more  convenient  to  use  but  present  difficulties,  1)  in  the 
attempt  to  determine  if  the  analytical  functions  satisfactorily  fit  the 
data  and  2)  in  determining  the  value  of  using  the  analytical  function 
beyond  the  range  of  the  data. 

In  the  effort  to  determine  the  probability  densities  of  the  HICAT  root -mean - 
square  gust  velocities  the  following  assumptions  were  made: 

a)  The  sample  size  was  adequate  for  only  the  root -mean-square  gust 
velocities  computed  from  the  spectra  with  the  truncation  limits 
130  ft.  and  2,000  ft. 

b)  The  analytic  function  is  given  by 

P  (x)  =  s  /  (2/tt)  exp  j  (2) 

where  x  represents_the  root-mean-square  gust  velocity  and  hence  can 
only  he  positive,  x  is  the  minimum  value  of  the  root-mean -square 
gust  velocity  for  the  air  to  be  considered  turbulent,  and  s  is  the 
standard  deviation  of  the  root-mean-square  gust  velocity. 

c)  x  =  0.5  ft/sec  for  the  vertical  component  and  x  =  0.5m^  +  9^  ft/sec 
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♦Numbers  in  parentheses  are  for  all  observations  given  equal  weight.  All  other  numbers  are  based,  upon 
weighting  observations  according  to  3ecgth  of  time  in  turbulence. 


for  the  lateral  component  and  x  =  0.5*g  +  0^  ft/sec  for  the  longi¬ 
tudinal  component.  The  values  of  bl,  ul,  6  and  are  taken  from 
Table  IV. 

d)  The  standard  deviations  of  the  root -mean -square  gust  velocities 
derived  from  the  HICAT  data  are  as  given  in  Table  V. 

•TABLE  V 


Standard  Deviations  of  RMS  (X  =  2,000  ft) 
Gust  Velocities 


Vertical 

Component 

Lateral  Longitudinal 
Component  Component 

Number 

Seconds 

Water 

(0.597)  0.456 

(0.801)  o.6o4 

(O.5O5)  0.409 

6210 

Flatland 

n 

(0.564)  0.608 

(O.634)  0.792 

(0.480)  0.535 

16845 

V  t 

Low  Mountains 

(0.879)  O.858 

(0.859)  1.00 

(0.817)  O.915 

8137 

High  Mountains 

(1.47)  1.25 

(1.52)  1.27 

(1.42)  1.27 

19388 

All  Cases 

0.938 

O.960 

0-910 

45,000-49,900  ft. 

0.6l0 

O.69O 

0.484 

5510 

50,000-54,900  ft. 

O.963 

1.027 

0-930 

20095 

55,000-59,900  ft. 

O.897 

0.830 

0.832 

15445 

>  60,000  ft. 

1.09 

1.20 

1.19 

9900 

*  Figures  in  parentheses  are  based  upon  using  equal  weights  for  all  turbulence 
observations.  The  other  figures  arc-  based  upon  calculations  weighting  each 
observation  by  the  time  In  turbulence. 


The  standard  deviations  based  upon  the  weighting  by  the  time  in  turbulence 
are  the  recommended  ones  because  weighting  each  run  equally  distorts  the 
contributions  of  the  short  runs. 

The  underlying  topography  appears  to  have  a  significant  effect  upon  the 
standard  deviation  of  the  root -meah-square  gust  velocity.  The  figures 
presented  In  Table  V  indicate  an  increase  in  the  standard  deviation  with 
an  increase  in  surface  roughness.  The  standard  deviation  for  the  lateral 
component  is  the  largest  for  all  of  the  categories  of  underlying  topography. 
For  the  water  and  flatland  categories  the  longitudinal  component  has  the 
smallest  standard  deviation  but  as  the  terrain  becomes  rougher  the  longl- 
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tudinal  component  Increases  relative  to  the  vertical  component. 

The  standard  deviations  for  all  three  components  are  smallest  in  the  45*000 
ft.  to  50,000  ft.  altitude  band  and  greatest  for  the  altitudes  above  60,000 
ft.  This  increase  in  the  standard  deviation  with  altitude  does  not  hold  in 
the  50*000  -  60,000  ft.  altitude  range.  Almost  all  turbulence  above  60,000 
ft.  altitude  occurred  over  mountains  and  vas  generally  more  intense  than  the 
average  turbulence  encounter  below  60,000  ft.  The  data  are  inadequate  to 
make  it  possible  to  definitely  separate  the  altitude  effect  from  the  terrain 
effect. 


The  probability  densities  of  the  root -mean -square  gust  velocities  were  com¬ 
puted  from  Equation  2  and  the  results  are  shown  in  Figures  1-5*  It  is 
important  to  note  that  these  probability  densities  relate  to  the  truncated 
spectra  and  not  to  the  total  spectra.  The  probability  densities  for  the 
root -mean-square  gust  velocities  derived  from  the  total  spectra  are  dis¬ 
cussed  in  Section  IV. 

Figures  6  through  11  consist  of  exceedance  curves  of  RMS  (X  =  2,000ft)  for 
each  of  the  three  components  and  for  the  four  categories  of  underlying 
terrain.  The  straight  lines  are  exceedance  curves  based  upon  the  standard 
deviation  computed  from  the  HICAT  data  and  the  assumption  that  the  distri¬ 
butions  are  normal.  The  standard  deviations  were  computed  by  using  both 
plus  and  minus  values  of  RM3  (X  =  2,000  ft)  -  x.  These  are  the  values  given 
in  Table  V.  The  data  points  fit  the  normal  distribution  reasonably  well 
except  for  the  cases  over  water  and  the  highest  values  of  RM3  (X  =  2,000  ft) 
for  the  mountain  cases.  There  are  insufficient  data  at  this  time  to  make 
possible  a  definitive  conclusion  relating  to  these  departures  from  normal. 
The  sample  size  is  relatively  small.  In  this  report  the  assumption  will  he 
made  th»t  the  distribution  of  HMS  (X  =  2,000  ft)  is  normal. 


If  the  distributions  of  KMB  gust  velocities  are  assumed  to  be  normal  for  the 
classifications  of  terrain  or  altitude  then  the  distribution  for  all  cases 
put  in  one  group  cannot  be  normally  distributed  because  the  sum  of  two  normal 
distributions  with  different  standard,  deviations  is  not  a  normal  distribution. 
The  probability  density  distribution  for  all  cases  then  would  he  represented 
by  « 

U*2 


P  (d)  =  P_ 


JTi 

1 TT  8, 


28, 


+  p„ 


JTi 

ITT  S_ 


2s„  +  - 


(3) 


where  the  P^,  Pg  -  -  -  -  P^  would  be  the  proportion  of  time  in  turbulence 
for  either  the  terrain  or  altitude  classifications  and  the  s^,  sg  -  -  -  s^ 
are  the  appropriate  standard  deviations. 


Effect  of  the  Definition  of.  Turbulence  Upon  the  Standard  Deviations  of  the 
Root -Mean -Square  Gust  Velocities. 

The  results  given  in  Table  V  were  based  upon  x  =  0.5  ft/sec  for  the  vertical 
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component  of  HMB  (X  =  2,000  ft).  If  tbe  definition  of  the  lower  bound  of 
turbulence  for  this  component  is  changed  to  x  =  1.0  ft/sec  then  the  standard 
deviations  of  the  weighted  values  are  as  given  in  Table  VI. 

r 

TABLE  VI 

Comparison  of  Standard  Deviations  of  RMB  (X  3  2,000  ft)' 

When  x  =  0.50  and  x  =  1.00  ft/sec  for  the  Various  Components 


Vertical  Component 

Lateral  Component 

Longitudinal  Component 

Underlying 

■pp 

X 

1  ; 

1 _ : 

X 

Terrain 

0.50  ft/sec 

Water 

0A56 

0.409 

O.60I* 

mm 

mm 

Flatland 

0.608 

1  0.41*5 

0.792 

mm 

Low 

Mountains 

ft 

CO 

• 

0 

0.1*82 

1.00 

■ 

■ 

0-534 

High 

Mountains 

• 

ro 

0.91*6 

1.27 

m 

O.969 

These  data  suggest  that  the  standard  deviation  showi  lees  variation  with 
different  x  values  than  with  classifications.  In  this  report  x  =  0.5  ft 
was  used  because  an  examination  of  the  data  indicated  that  this  value  appears 
to  be  in  closer  agreement  with  the  definition  of  turbulence  used  by  Steiner 
(5).  *  ; 
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igure  1.  Comparison  by  Component  and  Altitude  of  Probability  Density  of 
RMS  (X  =  2,000  ft)  -  x. 
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Figure  2.  Comparison  by  Altitude  and  Component  of  Probability  Density  of 
RMB  (X  »  2,000  ft)  -  x. 
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Figure  3.  Comparison  by  Component  and  Underlying  Topography  of  the  Probability 
Densities  of  »6  (\  =  2,000  ft)  -  x. 
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MM  MT  VELOCITY  (FT/SCQ 


Figure  4.  Comparison  by  Underlying  Topography  and  Component  of  the  Probability 
Densities  of  RMS  (X  -  2, 000  ft)  -  X. 
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Figure  5.  Comparison  by  Component  of  the  Probability  Densities  of  All 
Cases  of  RMS  (X  =  2,000  ft)  -  x. 
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Figure  6.  Probability  of  RMB  (X  =  2,000  ft)  -  x  Equalling  or  Exceeding 
Given  Magnitudes 
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Figure  7.  Probability  of  FMB  (X  -  2,000  ft)  -  x  Equalling  or  Exceeding 
Given  Magnitudes 


Figure  8.  Probability  of  RMB  (X  =  2,000  ft)  -  x  Equalling  or  Exceeding 
Given  Magnitudes 
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Figure  9.  Probability  of  RMB  (X  -  2,000  ft)  -  x  Equalling  or  Exceeding 
Given  Magnitude b 
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SECTION  HI 


SCALE  LEHGTES  OF  HIGH  AIT  ETUDE  CLEAR  AIR  TURBUIEKE 
)hthematlcal  Representation  of  Spectra 

Current  procedures  for  computing  power  spectral  densities  of  Measured  gust 
velocities  yield  numerical  values  that  define  irregularly  shaped  power 
spectral  density  curves.  Smooth  curves  obtained  by  "fairing"  or  by  "least 
squares  fit"  are  used  to  represent  the  principal  features  of  these  power 
spectral  density  curves.  The  mathematical  expressions  that  are  most  fre¬ 
quently  proposed  to  represent  these  smoothed  power  spectral  density  curves 
contain  three  parameters.  One  of  these  parameters,  the  root -mean-square 
gust  velocity,  9,1s  a  measure  of  the  Intensity  of  the  turbulence  and  has 
the  effect  of  moving  the  curve  upward  or  downward  In  a  direction  parallel 
to  the  ordinate.  The  second  parameter,  a,  determines  the  slope  of  the 
curve  at  the  high  frequency  end.  The  third  parameter,  designated  as  the 
scale  length,  L,  is  a  shape  parameter  that  defines  the  rate  of  change  of 
slope  at  the  low  frequency  portion  of  the  curve* 

Crooks,  Soblit,  and  Mitchell  (2)  discuss  In  detail  four  families  of 
functions  that  have  been  proposed  to  represent  power  spectral  density 
curves.  The  first  family  Includes  both  the  Van  Harman  spectrum  (m  =  -  5/3)  — 

and  the  Dryden  spectrum  {  a  =  -2)  as  special  cases.  The  general  form  of  the 
equation  for  this  family,  designated  the  Taylor-Bullen  family,  1b 


1(0) 


.Sli  [i  +  e  (  n  +  i  )  b2  (Op2  ] 

f  [1  .  b?  (CS,)2  ]“  *  V2 


where  n  and  b  are  constants. 


w 


The  second  family,  designated  the  Taylor-Bullen  longitudinal  family,  is 
defined  by 


#(n)  = 


2  or  L 


[1  +  b2  (n  L)2]n  +  2 

The  third  family,  designated  the  sharp-knee  family,  is  defined  by 


(5) 


t(n) 


Constant 
1  +  (fl  L)”m 


(6) 


and  the  fourth  family,  proposed  by  Hoblit  and  designated  the  mild -knee 
family,  is  defined  by 


#(n) 


Constant 

(1  +  n 


(7) 


♦ 


% 
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where  the  Constant  = 


i  t  2 

4  ir  L  cr 

3 

Figures  12  a-e  illustrate  the  curves  defined  by  Taylor-Bullen,  sharp-knee 
and  the  mild -knee  equations  form  =  -5/3  and  for  five  values  of  the  scale 
length,  L.  Inspection  of  these  curves  and  of  the  equations  leads  to  the 
following  conclusions  relating  to  the  interpretation  of  the  HICAT  spectra 
given  by  Crooks,  et  al,  (1,  2). 

a)  For  the  case  of  the  low  frequency  truncation  line  at 
2.5  x  lO-4  cycles/ft.,  there  ic  a  significant  difference 

,  between  the  equations  when  L  <  1,000  ft.  For  1,000  <  L 

<  4,000  ft.  the  "mild-knee"  equation  differs  signifi¬ 
cantly  from  the  other  equations.  For  L  >  4,000  ft.  the 
difference  between  any  two  of  the  equations  is  insigni¬ 
ficant. 

b)  For  the  low  frequency  truncation  line  at  10~^  cycles/ft., 
there  is  a  significant  difference  between  the  "mild -knee" 
equation  and  the  other  equations  for  L  <  8,000  it.  When 
L  >  8,000  ft.  there  is  no  significant  difference  between 
the  equations. 

Crooks  et  al.,  (1,2)  computed  the  root-mean-square,  gust  velocities  for  the  low 
frequency  truncation  line  at  10”3,  2.5  x  10"  ,  10“  ,  and  5  x  10 "5  cycles/l't. 
There  were  relatively  few  cases  for  the  truncation  line  at  5  x  10  cycles/ft. 

If  the  hypothesis  is  made  that  the  power  spectral  density  curves  are  adequately 
represented  by  any  of  the  equations  listed  above  and  that  the  slope  at  the  high 
frequency  (7.4  x  10”-*  cycles/ft. )  end  of  the  curve  is  known,  then  the  -value  of 
the  scale  length,  L,  can  be  related  to  the  ratio  of  RMS  (2.5  x  10”^)/RM5  (10"3 
cycles/ft.)  and  also  of  RM3  (10  cycles/ft.)/  R*6  (10“3  cycles/ft.).  Tables VII 
-VIII  and  Figures  13  and  14  give  the  values  of  these  ratios  for  various  values 
of  the  slope,  m,and  scale  length,  L.  A  detailed  analysis  of  the  method  used 
to  calculate  scale  lengths  from  RM9  ratios  is  presented  in  Appendix  I. 

Figures  13  and  14  and  Table  DC  indicate  that  for  values  of  -5/3 1  the 

difference  in  the  ratios  of  the  root -mean -square  gust  velocities  for  scale 
lengths  averaging  2,000  ft.  and  4,000  ft.,  using  the  mild-knee  equation,  is 
relatively  large.  Consequently,  if  the  mild-knee  equation  adequately  re¬ 
presents  the  HICAT  spectra,  then  only  gross  determinations  of  the  scale 
length  can  be  made. 

Scale  Lengths  for  the  Total  HICAT  Sample 

Scale  lengths  representative  of  the  total  HICAT  sample  may  be  obtained  by 
several  methods.  Crooks,  Hbblit,  and  Mitchell  (2)  obtained  composite  spectra 
by  two  different  methods  and  then  estimated  the  scale  lengths  by  visually 
comparing  the  composite  spectra  with  the  curves  defined  by  Equations  4-7* 

The  first  set  of  composite  spectra  were  obtained  by  first  preparing  cumulative 
probability  curves  foi  the  psd  values  at  each  of  various  frequencies.  Each 
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Figure  12  (Cent. ) 


value  was  weighted  according  to  the  length  of  the  turbulent  region.  Then 
probability  levels  were  selected  and  the  ped  values  corresponding  to  these 
probability  levels  were  plotted  as  a  function  of  frequency.  In  this  way 
composite  spectra  were  obtained  that  represented  the  0.5,  0.1  and  the  0.01 
probabilities,  respectively.  This  method  of  obtaining  composite  spectra  has 
the  advantage  of  clearly  Indicating  any  differences  that  might  exist  In  the 
spectra  representing  turbulence  of  different  Intensities.  The  second  set  of 
composite  spectra  was  defined  by  the  average  of  the  normalized  spectra. 

The  results  obtained  by  Croaks,  et  al.,  are  suomarlzed  in  Table  10. 

TABLE  Vn.  BMS  (2.5  x  ID"1*  cycles/ft.  )/HHS  (10"3  cycles/ft.)  for  the  Mild -Knee 
Equation  (Equation  j). 
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-1.6 

-1.67 

-1.7 

-1.8 

-1.9 

HESS 

500  ft. 
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1.459 
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1,000 

1.361 

1.407 
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1-704 

2,000 

1.404 

1.458 

1.584 

1.634 

1.657 

1.735 

1.821 

3,000 

1.421 
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1.544 

1.615 

1.669 

1.693 

1.778 

I.870 

1.970 

4,000 

1A31 

1.491 

1.5580 

1.632 

1.688 

1.713 

1.801 

1.897 

2.002 

5,000 

1.437 

1.499 

1.567 

1.642 

1.670 

1.725 

1.816 

1.915 

2.022 

6,000 

1.441 

1.504 

1-573 

1.650 

1.708 

1.734 

1.826 

1.927 

2.036 

8,000 

1.446 

1.510 

T.581 

1.659 

1.719 

1.745 

1.839 

1.942 

2.055 

1.450 

1.514 

1.586 

1.665 

1.725 

1.752 

1.848 

1.952 

2.066 

1.454 

1.520 

1-593 

1.673 

1.734 

1.762 

1.859 

1.965 

2.081 
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VIII.  bmB  (10  ^  eycles/ft)/  RMS  (10  ^  cyclea/ft)  for  the  Mild -Knee 
Equation  (Equation  7) 
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All  of  the  composite  spectra  obtained  from  either  the  cumulative  probability 
procedure  or  the.  averaging  procedure  with  the  spectra  truncation  line  taken 
at  fl  *  2.5  x  10”  cyclea/ft  (X  *  4,000  ft)  indicate  an  infinite  scale  length. 

The  scale  lengths  aay  also  be  determined  through  the  use  of  the  computations 
given  in  Tables  VII  and  VIII  and  Figures  13  and  14  combined  with  the  values  of  the 
root -mean-square  gust  velocities  published  by  Crooks,  et  al. ,  (l,  2).  If  the 
slopes  of  the  average  spectra  (Table  X)  are  used,  scale  lengths  of  1,200, 

6,000  and  3,000  ft.  for  the  vertical,  lateral,  and  longitudinal  components , 
respectively , are  indicated.  If  only  the  sets  of  spectra  are  used  for  which 
the  vertical  component  of  the  BIB  (5  x  10 cycles/ft)  is  equal  to  or  greater 
than  1.0  ft/sec,  then  the  slopes  are  found  to  be  as  indicated  for  the  0.5 
probability  level  and  the  scale  lengths  are  2,000  ft,  3,000  ft,  and  “  for  the 
vertical,  lateral  and  longitudinal  components, respectively. 

In  Figure  15  composite  spectra  and  their  associated  scale  lengths  and  slopes 
obtained  by  yet  another  method  are  presented  for  the  three  components  of  the 
gust  velocity  using  the  total  HKAT  sample.  The  composite  values  were  arrived 
at  by  (l)  averaging  the  logarithms  of  the  ordinate  values  at  various  fre¬ 
quencies  and  (2)  taking  the  antilog  of  the  averages.  The  advantage  of  this 
method  over  one  which  uses  direct  averages  of  the  ordinates  is  that  the  latter, 
because  of  distortion  due  to  thexlogarithm  scale,  tends  to  give  unequal  velgat 
to  spectra  with  less  bending  at  low  frequencies  whereas  composite  spectra  de¬ 
rived  from  averaging  logarithms  more  closely  resemble  the  median  value. 

8cale  lengths  of  over  8,000  feet  are  obtained  by  this  method  of  averaging  the 
logarithms  of  the  ordinates.  Spectra  with  root-mean -square  gust  velocities 
(truncated  at  5  *  10  cydes/f+)  >  1.0  ft/sec  and  extending  to  10  *  cycles/ 
ft  were  used  in  obtaining  the  composite  spectra. 

Figure  l6  illustrates  the  results  obtained  when  the  data  used  to  determine  the 
spectra  for  the  vertical  component  given  in  Figure  15  are  divided  into  two  in¬ 
tensity  classes.  The  composite  spectrum  representing  the  high  intensity  cases 
has  a  steeper  slope  and  a  shorter  scale  length  than  the  spectrum  for  the  low 
intensity  cases  but  the  difference  appears  insignificant. 
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Figure  13,  Scale  Lengtn,  L,  as  a  Function  of  the  Ratio  of  Root -Mean -Square 


•  Gust  Velocities  and  Slope ,m,  for  the  Mild -Knee  Equation.  (Shaded 

*  Areas  Indicate  Effects  of  2#  and  5$  Uncertainties  in  the  Abscissa 

|  and  Ordinate  Values) 
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Figure JA.  Scale  Length,  L,  as  a  Function  of  the  Ratio  of  Root -Me  an -Square 
•  !  Gust  Velocities  and  Slope., m.  for  the  MLld-Khee  Equation.  (Shaded 
Areas  Indicate  Effects  of  2>  and  5^  Uncertainties  in  the  Abscissa 
and  Ordinate  Values) 


TABLE  IX.  Range  In  Computed  Value*  of  Scale  length  for  Varloua  Errors  in 
a  and  Ratio  of  the  Root -Mean -8quare  Oust  Velocities. 


Error  in  RMS 
(2.5-x  10“4)/BM5 

do"3) 

Range  in  L 
(Mean  t  «  2,000  ft) 

Range  in  L 
(Mean  L  *  *,000  ft) 

..±2* 

1700  to  2600  ft 

2800  to  6000  ft 

±5* 

9QO  to  6000  ft 

1500  to  * 

Error  in  m 
(Mean  a  =  -I.67) 

±  2* 

1700  to  2600  ft 

3000  to  8000  ft 

±5* 

11*00  to  1*500  ft 

2000  to  * 

Error  In  Ratio  of 

RM3  Values  Combined 
with  error  in  m 

±  2* 

1500  to  35OO  ft 

2200  to  15000  ft 

i  5* 

750  to  * 

1000  to  • 

Error  in  RMS  _ 

(10~  ;/RM3  (10~3) 

+  1 

1700  to  2500  ft 

3300  to  5000  ft 

±5* 

ll*00  to  3000  ft 

2500  to  9000  ft 

Error  in  m 
(Mean  m  =  -1.67) 

* 

i  2^ 

1700  to  21*00  ft 

3200  to  5600  ft 

±  'll 

1500  to  3500  ft 

2500  to  13000  ft 

Error  in  Ratio  of 

RMS  Values  Combined 
with  error  in  tn 

1500  to  3000  ft 

noo  to  6000  ft 

2600  to  7000  ft 

1700  to  • 
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INVERSE  WAVELENGTH  -  <c*» ) 

Figure  13.  Spectra,  Slopes  and  Scale  Lengths  for  the  Total  HICAT  Sample. 


INVERSE  WAVELENGTH  -  ( c*f ) 

Figure  16.  Comparison  of  Normalized  Composite  Spectra  Separated  by  Intensity 
of  the  Turbulence. 
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Scale  Lengths  of  Composite  Spectra  Representing  Topographic  and  Altitude  i 
Categories.  * 

Three  netbods  were  used  to  deterw&ne  the  scale  lengths  of  the  cospoeite  -j 
spectra  representing  the  topographic  and,-  altitude  categories.  The  results; 
are  presented  In  Table  XI.  In  the  first  sethod  the  composite  spectra  verer  i 
obtained  by  averaging  the  logarithms  of  the  individual  spectra  ordinate  •? 
values  as  explained  earlier  in  the  section.  These  spectra  are  illustrated} 
in  Figures  17  and  18.  In  the  second  and  third  netheds  scale  lengths  sere  ^ 
obtained  toy  comparing  calculated  values  of  EMB  (2.5  *  10  cycles/ft  )/HC  1 
(lO'J'  cycles/ft)  and  KMB  (lO*4  cycles/ft)/  EMB  (10"3  cyclae/ft)  with  the  1 
values  for  the  aild-knee  equation  listed  in  Tables  VH  an3  VH2.  A  detailed 
description  of  these  aetbods  is  found  in  Appendix  !•  j 


TABIZ  XI.  Beale  Lengths  of  Coapoelte  Spectra  ; 

U 


Category 

Vertical  Component 

Lateral  Component 

Figure 

17  A  18 

IjSffl 

ca 

Table 

vin 

Figure 
17  4  18 

Table 

VH 

Table 

vni 

Figure 
17  &  l8 

Table 

vni 

Altitude 
<  59,000  Ft 

8,000  ft 

a» 

CM 

16,000 

8,000 

CO 

16,000 

1 

B 

Altitude 
>  59,000  Ft 

8,000 

• 

B 

16,000 

8,000 

mb 

16,000 

B 

B 

Flatland 

4,000 

a 

4,000 

8,000 

■1 

12,000 

40,000( 

MB 

HI 

Low 

Mountains 

16,000 

■ 

00 

4,000 

3,500 

5,500 

Wm 

OB 

MB 

High 

Mountains 

4,000  ' 

m 

15,000 

8,000 

MM 

8,000 

16,000 

/ 

mb 

MB 

Computations  of  the  scale  lengths  can  be  sede  using  the  correlation  function 
sethod  described  by  Houbolt  (5).  This  sethod  requires  data  which  is  essentially 
unfiltered.  The  data  sasple  available  from  HEAT  records,  however,  was  bigh- 
Pmi  filtered  and  d»racterised  by  rapidly  decaying  correlation  functions  re¬ 
sulting  in  scale  lengths  about  one  order  of  nagnitude  lower  than  what  was  de¬ 
rived  using  other  techniques. 

In  the  following  section  the  scale  lengths  given  by  the  composite  spectra 
(Figures  15  to  10)  will  be  used  in  converting  the  root -ae an -square  gust 
velocities  obtained  fres  the  truncated  spectra  to  those  useful  in  aircraft 
gust  load  calculations* 
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Figures  l8a,b,c»  Caapoeite  Spectra  for  Altitude  Categories 


SECTION  IV 


STANDARD  DEVIATION  OF  THE  ROOT -MEM -SQUARE  GUST  VELOCITY 


Introduction 


Alrcrf ft  gust  load  calculations  are  based  in  part  upon  the  probability  density 
distribution  of  the  root-mean-square  gust  velocity.  This  probability  density 
distribution  i*  conventionally  define^  as  the  sum  of  two  normal  density  distri¬ 
butions  weighted  by  P^  and  P^, 

p  fa)  =  ^ exp  (fi/  2bi2)  + 

(  exp  2  b22)  (8) 

where  P,  and  P?  are  the  proportions  of  total  flight  time  in  "non-storm"  and 
"storm"-1  turbulence  respectively;  b^  and  bg  are  the  standard  deviations  of  the 
root -mean-square  gust  velocity^  a~,  in  non-storm  turbulence  and  storm  turbu¬ 
lence  respectively.  am  m  a  -  x;  °  is  the  true  root-mean-square  gust  velocity 
and  x  is  the  minimum  value  of  17  to  be  classified  as  turbulence.  The  root-mean- 
square  gust  velocity  is  equal  to  the  square  root  of  the  area  under  the  complete 
power  spectral  density  curve.  Power  spectral  density  curves  computed  from 
measurements  have  truncation  limits  that  are  a  function  of  the  length  of  the 
turbulence  record,  instrument  capabilities  and  the  computation  technique.  If 
the  power  spectral  density  curves  computed  from  the  gust  velocity  measurements 
are  found  to  be  adequately  represented  within  the  frequency  range  of  the 
measurements,  by  an  equation  (e.g.  one  of  Equations  4  llucugii  i  /  ujuen  uuxu 
equation  may  be  used  to  determine  the  area  under  the  curve  with  the  frequency 
•limits  set  at  0  and  «,  providing  an  assumption  ib  made  that  the  equation 
adequately  represents  the  spectra  at  frequencies  outside  the  range  of  the 
measurements. 


This  assumption  suggests  the  following  example.  Consider,  two  observations 
of  turbulence  with  equal  root-mean-square  gust  velocities  computed  from 
spectra  truncated  at  5  -  10  ,  cycles/ft.  (X  =  2,000  ft).  In  one  case  the 
turbulence  region  is  Just  long  enough  to  compute  the  spectrum  out  to  this 
truncation  limit.  In  the  second  case  the  turbulence  region  is  long  enough 
to  provide  a  spectrum  to  a  truncation  limit  of  10 cycles/ft.  (X  =  10,000 
ft),  fJuppose  that  the  two  spectra  have  the  same  slope  and  curvature  up  to 
5  x  104  cycles/ft.  (X  -  2,000  ft)  and  hence  can  be  represented  by  the 
equation  within  the  common  frequency  range.  In  such  a  circumstance 
the  physical  situation  in  the  turbulent  region  from  which  the  shorter 
record  was  obtained  is  assumed  identical  to  that  of  the  longer  record  * 
which  is,  of  course,  not  necessarily  so.  However,  for  purposes  of  this 
section  the  assumption  is  considered  valid. 


♦  f  - 

Estimating  b^'a  *  , 

In  Section  II  the  probability  density  distributions  were  determined  for  the 
measured  root-aean-square  gust  velocities  obtained  from  the  spectra  between 
the  truncation  limits  7.4  x  10 and  5  x  cycles/ft.  These  probability 
density  distributions  were  converted  into  probability  density  distributions 
of  through  use  of  the  following  proceedure: 

a)  The  "mild -knee"  equation  (Equation  7>  Section  IH)  was  assumed  to 
best  represent  all  the  power  spectral  density  curves.  It  should  be 
noted,  however,  that  for  L's  on  the  order  of  those  derived  from  the  . 
HICAT  spectra  (8,000  -  l6,000  ft)  the  choice  of  equations  has  only 
a  small  effect  on  the  results. 


b)  The  values  of  scale  length,  L,  and  slope,  m,  given  for  the  various 
composite  spectra  discussed  in  Section  111  were  applied  to  the  in¬ 
dividual  spectra.  Consequently,  the  individual  spectra  were 
assumed  to  differ  from  eac^  other  only  in  the  "scale  factor"  that 
indicates  the  relative  intensity  of  the  turbulence. 

• 

c)  Table  XII  was  used  to  convert  the  root -me  an -square  gust  velocities 
obtained  from  the  truncated  spectra  to  the  root -mean-square  gust 
velocities  obtained  from  Equation  7  between  the  limits  0  and  08 
cycles /ft. 

*  .  . 

The  Justification  for  the  assumption  made  in  (a)  was  that  those  spectra  Ob¬ 
tained  from  the  long  turbulent  regions  were  best  represented  by  the  "ndld- 
knee"  equation.  As  the  discussion  in  Section  III  indicated,  no  definite 
choice  between  Equations  4  through  7  xould  be  made  for  the  spectra  obtained 
from  the  relatively  short  turbulent  regions. 


The  Justification  of  the  assumption  made  in  (b)  was  that  the  scale  of  ..tur¬ 
bulence.  L,  could  not  be  determined  for  the  spectra  with  the  truncation  at 
5  x  10“^  cycles/ft.  and  hence  values  for*  L  had  to  be  assumed.  The  L's  gi'ucr; 
by  the  composite  spectra 'appeared  to  be  the  most  reasonable.  The  slopes  of 
the  composite  spectra  were  used  so  as  to  he  consistent  with  the  scale 
lengths. 


The  standard  errors  of  the  b  's  listed  in  Table  XIII  are  estimated  to  be  at\ 
least  of  the  values  listed.  The  uncertainties  in  the  values  of  the 
b.'s  arise  principally  because  power  spectral  density  curves  computed  from 
the  gust  velocity  measurements  cover  a  relatively  limited  frequency  range. 
Errors  on  the  evaluation  of  the  b1's  arise  because 

a)  The  values  of  s,  (standard  deviation  of  BMB  (X  ■  2,000  ft)  may  vary 
as  much  as  50^  as  the  definition  of  the  lower  boundary  of  turbulence 
is  changed  over  reasonable  limits.  In  the  preparation  of  Table  XIII 
turbulence  was  defined  to  exist  if  Uy  (5  x  10”  cycles/ft)  >  0.5  ft/ 
sec. 


a 
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TABLE  XU 


4  Ratio  of  a_  to  Root -Mean-Squar?  Giiat  Velocity 

1  i 

From  Spectra  Truncated  at  5  x  10 "  Cyclea/Ft  (Based  on  Mild -Knee  Equation) 

X 

\ 


L 

■ 

m 

B 

D 

I.67 

/ 

B 

1.8 

B 

m 

500  Ft 

1.29  - 

1.33 

1.38 

1.43 

1.46 

1.48 

1.54 

1.60 

1 .66 

1,000 

1.42 

1.49 

1-57 

1.66 

1.73 

1.76 

1.86 

1.98 

2.11 

Mg 

1-59 

1.70 

1.83  - 

1.98 

2.10 

2.15 

2.34 

2.55 

2.80 

3,000 

1.70 

1.85 

2.02 

2.21 

■1 

2.44 

2.70 

m 

3.35 

1.  AAA 

*»  nQ. 

JLm  |W 

IS 

O  lirt 

a-*  -tv/  1 

O  CO 
•-m  J7 

2.68 

3-07 

■3 

- 

5,000 

1.85 

2.05 

2.28 

2.56 

2.78 

2.89 

3.27 

3.71 

'  4.23 

6,000 

1.91 

2.13 

2.39 

2.70 

2.95 

3«oi 

3.50 

4.01 

4.62 

8,000 

2.00 

2.26 

2.57 

2.94 

3.24 

3.38 

3.92 

4.55 

5.30 

10,000 

2.08 

2.37 

2.72 

3.14 

B 

3.65 

4.27 

5.01 

5.90 

15,000 

2.23 

2.58 

3.00 

3.54 

3-99 

4.20 

5.00 

5-99 

B 

\ 


b)  The  scale  length,  L,  can  only  be  crudely  estimated  for  the  spectra 
obtained  f ran  short  records.  Here  it  was  assumed  that  these 
spectra  have  the  same  distribution  of  L's  as  the  spectra  from  long 
records. 

c)  Various  methods  of  determining  composite  spectra  yield  different 
values  of  L  and  m.  The  values  given  in  Table  XIII  are  considered 
a  best  estimate. 

If  the  definition  of  the  lover  limit  of  turbulence  that  is  given  in  (a)  is 
accepted  then  the  lover  limit  of  for  turbulence  becomes  a  function  of  L 
and  m.  These  lover  limits,  x  (see  equation  8), are  also  listed  in  Table  XIII. 
In  Section  V  the  lover  limit  of  turbulence  used  to  define  the  ratio  of  tur¬ 
bulent  flight  miles  to  total  flight  miles  is  defined  in  terms  of  the  eg 
normal  acceleration,  vhere  frequently  occurring  peaks  of  +  0.10g  ere  equated 
to  a  KH3  (X  =  2,000  ft)  of  0.5  ft/sec. 
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SECTION  V 


RATIO  OF  TURBULENT  FLIGHT  MILEfS  TO  TOTAL  FLIGHT  MILES 
Definition  of  Turbulence 


When  an  aircraft  that  la  flying  through  a  portion  of  the  atmosphere  in  which 
there  are  nc  clouds  undergoes  accelerations  that  cannot  be  directly  attri¬ 
butable  to  the  movement  or  setting  of  the  control  surfaces,  the  aircraft  is 
said  to  be  in  clear  air  turbulence.  The  accelerations  ot  the  aircraft  are 
functions  of  the  weight  and  characteristics  of  the  aircraft  in  addition  to 
the  atmospheric  gusts.  If  the  ratio  of  the  turbulent  flight  miles  to  the 
total  flight  miles  is  to  be  computed,  definitions  of  turbulent  and  total 
flight  miles  are  required.  -The  definition  of  the  turbulent  flight  miles 
oust  include  statements  that  give  (l),  the  lower  limit  of  accelerations 
that  are  considered  to  be  turbulence,  (2),  the  frequency  interval  of  interest, 
and  (3),  a  quantitative  evaluation  of  the  duration  of  the  turbulence.  If  the 
turbulent  regions  are  not  randomly  distributed  in  space  and  time,  the  ratio 
of  turbulent  flight  miles  to  total  flight  miles  is  also  a  function  of,  (l), 
the  distribution  and  total  number  of  flight  miles  by  altitude,  season,  and 
topographic  region,  (2),  turbulence  search  or  avoidance  procedures.  If 
successful,  and  (3),  pattern  flying  through  a  known  turbulent  region. 

In  Section  II,  for  the  purpose  of  determining  the  probability  density  distri¬ 
bution  of  the  vertical  component  root -mean -sauare  gust  velocity,  turbulence 
was  defined  to  exist  when  RMS  (X  =  2,000  ft)  was  equal  to  or 
greater  than  0.5  ft/sec.  TBis  definition  of  clear  air  turbulence  is  in 
apparent  conflict  with  the  definition  given  in  terms  of  the  aircraft  acceler¬ 
ation.  For  the  purpose  of  determining  the  ratio  of  the  turbulent  flight 
miles  to  the  total  flight  mile's,  the  definition  of  turbulence  given  below  in 
terms  of  acceleration  of  the  aircraft  was  used  because  the  gust  data  were  in¬ 
complete  from  the  HICAT  flight  program  and  non-existent  from  the  NASA  U-2 
flight  program.  In  addition,  the  turbulent  region  had  to  be  approximately 
8  miles  long  to  provide  adequate  data  for  a  spectra  of  minimum  usable  length. 
Thus,  if  the  RMS  gust  velocity  was  used  to  define  turbulence,  all  turbulent 
regions  shorter  than  8  miles  would  be  ignored. 

In  his  review  of  the  VGH  data  from  768,000  miles  in  the  altitude  range  40,000 
to  70,000  ft.  of  NASA  U-2  flights  Steiner  (  5  )  defined  turbulence  to 
exist  "whenever  the  accelerometer  trace  was  disturbed  and  contained  gust 
velocities  (presumably  "derived"  gust  velocities,  ”de>  greater  than  2  ft/sec." 
Crooks  (l)  defined  turbulence  to  exist  whenever  the  eg  normal  acceleration 
trace  showed  frequently  occurring  peaks  of  +  0.10  g.  The  beginnings  and 
endings  of  turbulence  were  given  to  the  nearest  5  seconds  and  turbulence  of 
less  than  ten  seconds  duration  was  counted  as  no  turbulence.  A  comparison 
of  the  mum  U.  values  of  the  HICAT  and  the  NASA  flights  indicates  that 
Crooks'  and  Steiner's  criteria  for  turbulence  are  approximately  equivalent. 

The  classification  of  the  HICAT  turbulence  and  non-turbulence  flight  miles 
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by  altitude,  underlying  topography  and  by  season  waa  done  subjectively. 

The  principal  difficulty  waa  associated  with  patterns  or  repeated  flights 
through  a  region  that  was  found  to  be  turbulent  pn  the  flight  through  the 
region.  In  general,  the  turbulence  was  found  to  vary  in  intensity  and 
duration  for  each  flight  through  the  region.  Approximately  3^000  flight  » 
miles  of  the  more  than  300,000  flight  miles  were  used  in  the  determination 
of  the  ratio  of  turbulent  flight  miles  to  total  flight  miles.  The  155*000 
miles  not  used  were  associated  with  pattern  flights,  test  flights  and  in¬ 
strument  failure. 

Ratio  of  Turbulent  Flight  Miles  to  Total  Flight  Miles  for  Categories  of 
Altitude,  Topography  and  Season. 

Table  XIV  lists  the  estimated  ratio  for  the  HICAT  flights  of  the  turbulent 
flight  miles  to  the  total  flight  miles  for  categories  of  altitude,  topo¬ 
graphy  and  season.  The  topographic  categories  were  defined  as  follows: 

a)  Flatland  -  relief  differences  less  than  3*000  ft 

b)  Low  Mountains  -  relief  differences  3*000  to  7*000  ft 

c)  High  Mountains  -  relief  differences  over  7*000  ft. 

Turbulence  was  presumed'  to  exist  if  the  eg  acceleration  trace  frequently  j 
exceeded  +  0.10te.  Table  XV  compares  the  ratios  given  by  Steiner  for  the 
NASA  U-2  flights,  the  HICAT  -  program  and  by  MIL-A-8861A  for  altilateJLntei-vals. 
Table  XVI  lists  the  ratios  of  the  turbulent  flight  miles  tb  totsjfiNfiht 
miles  that  are  recommended  by  the  authors  for  use  in  the  deslgnflrcwperation 
of  advanced  aircraft.  The  ratios  obtained  from  an  analysis  oi'^gJfcCAT  data 
were  Judged  to  be  high  mainly  because  of  the  search  procedures  used.  The 
NASA  flights  were  probably  distributed  more  nearly  to  the  expected  flight 
operations.  This  also  Justified  lowerirg  the  ratios  obtained  from  HICAT  data. 

The  proportion  of  the  total  flight  miles  that  are  turbulent  varies  signifi¬ 
cantly  with  the  definition  of  turbulence.  The  relative  magnitude  of  this 
variation  is  illustrated  in  Table  XVII.  This  table  alBO  clearly  indicates 
that  there  is  little  change  of  the  relative  frequency  of  occurrence  of 
turbulence  with  altitude  over  mountains  and  a  relatively  large  change  with 
altitude  over  flatland  and  water. 

The  numbers  of  HICAT  flight  miles  by  season  and  topography  are  shown  in 
Table  XVIII.  Flights  over  mountains  were  predominantly  made  in  the  winter. 
Flight  miles  were  approximately  equally  distributed  by  seas  chi  over  flat- 
land  and  water.  Table  XIX  indicates  the  proportion  of  the  total  flight 
miles  that  were  turbulent  in  categories  of  season  and  topography.  There 
was  relatively  little  turbulence  during  summer  and  autum  over  water  and 
flatland. 


TABLE  XIV 


Ratio  of  Turbulent  Flight  Niles  to  Total  Flight  Niles 


Category  " 

Total  Flight  Niles 

Ratio  Turbulent 

X 

65,000  -  69,900  ft 

l6,000  miles 

0.019 

60,000  -  64,900 

87,000 

0.014 

55,000  -  59,900 

107,000 

0.028 

50,000  -  54,900 

105,000 

0.047 

45,00a  -  49,900 

30,000 
*  • 

0.042 

45,000  -  69,900 

345,000 

0.031 

Water 

124,000 

0.027 

Flatland 

132,000 

0.027 

Lav  Mountains 

49,000 

0.036 

High  Mountains 

40,000 

0.049 

Winter  (Dec.  -  Feb.) 

108,000 

0.036 

Spring  (Mao*.  -  May  ) 

81,000 

0.039 

Sumer  (Jun.  -  Aug.) 

71,000 

0.019 

Autism  (Sept.  -  Nov.) 

85,000 

0.024 

42 
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TABLE  XV 


Comparison  of  Steiner  (NASA  U-2),  HEAT  and  MIL-A-8861A  Data 


Altitude 
*  Interval 

Steiner  (  ) 

- - 

HEAT 

MIL-A-8861A 

(1968) 

Ratio 

Turbulent 

Flight 

Miles 

Ratio  ' 
Turbulent 

Flight 

Miles 

Ratio 

Turbulent 

60,000-70,000  Ft 

576,000 

0.006 

103,000 

0.015 

OoOOl 

50,000-60,000  Ft 

141,700 

0.020 

212,000 

0.035 

0.002 

40,000-50,000  Ft 

49,500 

0.025 

30,000 

0.042 

0.005 

40,000-70,000  Ft 

768,100 

0.009 

345,000 

0.031 

TABLE  XVI 

Recommended  Ratios  of  Turbulent  Flight  Miles  to  Total  Flight  Miles 


Category- 


Ratio  Turbulent 


65,000-70,000  Ft 
60,000-65,000 


0.010 


55,000-60,000 


0.010 


0.020 


50,000-55,000 


45,000-50,000 

45,000-70,000 


0.030 


0.030 


0.020 


Water 


0.018 


Flatland 


0.018 


Low  Mountains 


High  Mountains 


0.022 


0.030 


Winter 


0.024 

0.024 


Spring 


o.oi4 

0.018 


Summer 


Autumn 


TABLE  XVII 


luaber  of  Occurrences  by  Terrain  of  Turbulence  with  Vertical  RMS 
Greater  than  Listed  Values  for  Given  Altitude  Bands 


Hater, 

Flatland  (100  Sables) 
Altitude  Band 

Mountain*  (120  SsatO.es) 
Altitude  Betrd 

Uv  RM5  (\  -  2,000  ft) 

<55K 

55-60K 

s60K  ft 

<55K 

55-&K 

*60Kft 

68 

26 

4 

32 

35 

33 

WKZEEE3EKM 

16 

6 

0 

20 

20 

MEM 

>1.5  ft/sec 

4 

4 

0 

13 

10 

2 

X 

0 

10 

5 

7 

PB289HI 

1 

0 

0 

7 

2 

2 

TABLE  XVIH 


Flight  Miles  by  Season  for  Four  Terrain  Types  (HICAT  Flights  54  to  265  ) 


Water 

Vinter 

Spring 

— 1 

Fall 

All 

Seasons 

PU.cra 

—  •# 

oc  Pftn 
J/l^ 

gglgj 

2.5,390 

Flatland 

28,740 

lew 

132,180 

Low  Mountains 

31,330 

5,580 

4,870 

7,240 

49,020 

High  Mountains 

22,390 

PPPpP 

6,370 

5,280 

39,640 

All  Terrain 

107,410 

81,920 

70,610 

84,790 

3^4,730 

TABLE  XIX  0 

Percentage  Flight  Miles  with  Turbulence  >  Light  (eg  peaka  >  + O.lOg)  and  Moderate  In 
Parenthesis  (eg  peaks  >  ♦  0.25g)  by  Season  for  Four  Terrain  Types  (HICAT  Flights  54  to  285) 


,  Hater 

Winter 

Spring 

Suaer 

Fall 

All 

Seasons 

jjjgai 

Q2SUI 

mogi 

BBSBI 

217(0.5) 

Flatland 

3*6(0. 9) 

g 

2*  7(0*5) 

Low  Mountains 

2. 9(0.4) 

3*5(o.6) 

3*6(l*7) 

High  Mountains 

KSISSUI 

IffiSSI 

HttStl 

4. 9(1*8) 

All  Terrain 

Biiffil 

EBSBI 

l^g| 

3.1(08) 

S3CTICH  VI 


TURBULENCE  ABOVE  THUNDERSTORMS 

Clear  Air  Turbulence  Above  Thunderstorms  in  the  5ICAT  Program 

In  the  HEAT  flight  program  nearly  one  tenth  of  all  the  clear  air  turbulence 
observed  during  the  periods  that  the  aircraft  was  flying  over  water  or  flat- 
land  areas  was  associated  with  flights  above  tnunderstorms.  During  the 
flight  program  the  pilots  definitely  identified  42  cases  of  flights  above 
thunderstorms.  Turbulence  was  observed  on  41  of  these  occasions.  Table  XX 
lists  some  of  the  characteristics  of  the  turbulence. 


TABLE  XX 

Characteristics  of  Turbulence  Above  Thunderstorms  , 

W 


Estimated  Altitude 
Above  Thunderstorm 

Number 
of  Cases 

Mean  Length 
.  of  Turbulent 
Region 

Ude 

(Mean  Max.) 

Mean  RMS 

Ude 

1,000  -  3,000  ft 

16 

22  nm 

7-0  ft/sec 

1.6o  ft/sec 

4,000  -  5,000  ft 

13 

17 

4.9 

A  AA 

u.  yc. 

6,000  -  9,000  ft 

12 

10 

2.5 

0.76 

>  9,000  ft 

1 

0  V  , 

- 

✓ 

The  mean  length  of  the  turbulent  region  for  those  cases  of  flights  1,000  to 
3,000  ft.  above  thunderstorms  Is  approximately  3 0J&  larger  than  the  estimates 
found  In  the  literature  of  the  mean  dlametdr  of  thunderstorms.  This  agree¬ 
ment  Is  considered  to  be  relatively  good  when  the  variability  of  the  size  of 
thunderstorm  and  the  small  HEAT  sample  are  considered.  Data  relating  to 
the  true  gust  velocities  are  not  given  lh  Table  XX  because  true  gust  velo¬ 
cities  were  only  given  in  four  cases.  In  general,  the  turbulent  regions  were 
net  long  pnough  to  provide  the  data  needed  for  the  power  spectral  density 
computations  over  a  useful  frequency  range.  ^The  decrease  in  the  length  of 
the  turbulent  regions  and  in  the  intensity  of  the  turbulence  as  Indicated  by 
the  derived  gust  velocity,  U,  ,  vith  Increasing  altitude  above  the  thunder¬ 
storm  appears  reasonable  although  the  results  may  be  fortuitous  because  of 
the  small  sample. 
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tabu:  xvii 


lumber  of  Occurrence*  by  Terrain  of  Turbulence  with  Vertical  RMS 
Greater  than  Listed  Values  for  Given  Altitude  Bands 


Water, 

Flatland  (100  SaqO.es) 
Altitude  Band 

Mountains  (120  SasgO.es) 
Altitude  Band 

Uv  RMS  (\  r  2,000  ft) 

<55K 

55-60K 

>60K  ft 

<55K 

55-60K 

>60Kft 

”*0.5  ft/sec 

68 

28 

4 

32 

35 

33 

>1.0  ft/sec 

16 

6 

0 

20 

20 

25 

>1.5  ft/sec 

4 

4 

0 

13 

10 

15 

>2.0  ft/sec 

2 

1 

0 

10 

.  5 

7 

>2.5  ft/sec 

1 

0 

0 

7 

2 

2 

TABLE  XVIII 


Flight  Miles  by  Season  for  Four  Terrain  Types  (HICAT  Flights  54  to  285  ) 


Water 

Winter 

Spring 

Summer 

Fall 

All 

Seasons 

24,950 

35,680 

37,670 

25,390 

123,890 

Flatland 

28,740 

34,860 

21,700 

46,880 

132,180 

Low  Mountains 

31,330 

5,580 

4,870 

7,240 

49,020 

High  Mountains 

22,390 

5,600 

6,370 

5,280 

39,640 

All  Terrain 

107,410 

81,920 

70,610 

84,790 

344,730 

TABLE  XIX 

Percentage  Flight  Miles  vith  Turbulence  >  Light  (eg  peaks  >  +  0.10g)  and  Moderate  in 
Parenthesis  (eg  peaks  >  +  0.25g)  by  Season  for  Four  Terrain  Types  (HICAT  Flights  54  to  285) 

_ r  _ _  _  _  _ 


Water 

Winter 

Spring 

Sumner 

Fall 

All 

Seasons 

KSiffil 

QggJJJI 

KKSSI 

KBSSM 

2.7(0. 5) 

Flatland 

KSEfiil 

JEitiEM 

KJSSSI 

^gggi 

2. 7(6. 5) 

Low  Mountains 

■B3QP1 

EE1SBSI 

EBilSBil 

3.6(i.7) 

High  Mountains 

K3332I 

6.6(1. 5) 

1 

4. 9(1.8) 

All  Terrain 

ESilSSI 

E22S8I 

1*9(0. 3) 

mmmm, 

3.l(o.S) 
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1 ,  SECTION  VI 

\ 

TURBULENCE  ABOVE  THUNDERSTORMS 
i 

•  Clear  Air  Turbulence  Above  Thunderstorms  In  the  RICAT  Program 

« 

In  the  HICAT  flight  program  nearly  one  tenth  of  all  the  clear  air  turbulence 
observed  during  the  periods  that  the  aircraft  vas  flying  over  water  or  flat- 
land  areas  vas  associated  with  flights  above  tmmderetonns.  During  the 
I  flight  program  the  pilots  definitely  identified  42  cases  of  flights  above 

;  thunderstorms.  Turbulence  vas  observed  on  4l  of  these  occasions.  Table  XX 

!  lists  some  of  the  characteristics  of  the  turbulence. 

TABLE  XX 

Characteristics  of  Turbulence  Above  Thunderstorms 


Estimated  Altitude 
Above  Thunderstorm 

Number 
of  Casea 

Mean  Length 
of  Turbulent 
Region 

Ude 

(Mean  Max.) 

Mean  RMS 

Ude 

1,000  -  3,000  ft 

l6 

22  nm 

7.0  ft/sec 

1.60  ft/sec 

4,000  -  5,000  ft 

» 

13 

17 

4.9 

0.92 

6,000  -  9,000  ft 

12 

10 

2.5 

O.76 

.  >  9,000  ft 

1 

0 

- 

- 

The  mean  length  of  the  turbulent  region  for  those  cases  of  flights  1,000  to 
3,000  ft.  above  thunderstorms  is  approximately  30 £  larger  than  the  estimates 
found  In  the  literature  of  the  mean  diameter  of  thunderstorms.  This  agree¬ 
ment  is  considered  to  be  relatively  good  vhen  the  variability  of  the  size  of 
thunderstorms  and  the  small  HICAI  sample  are  considered.  Data  relating  to 
the  true  gust  velocities  are  not’  given  in  Table  XX  because  true  gust  velo¬ 
cities  vere  only  given  In  four  cases.  In  general,  the  turbulent  regions  vere 
not  long  enough  to  provide  the  data  needed  for  the  power  spectral  density 
computations  over  a  useful  frequency  range.  The  decrease  in  the  length  of 
the  turbulent  regions  and  In  the  intensity  of  the  turbulence  as  Indicated  by 
the  derived  gust  velocity,  U.  ,  with  increasing  altitude  above  the  thunder¬ 
storm  appears  reasonable  altnSugh  the  results  may  be  fortuitous  because  of 
the  small  sample. 
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The  relative  frequei«v  *>''  T’-lgfil?  th.^ve  lUandarsio  tb  Lr.  the  '•'ju.'.YT  progrs,": 
may  oe  high  because  of  _h-.  turbvier.et.  -.tee  slog  program  and  because 

the  bases  and  lice  of  operation  were  tfeliterateii  •  :>v,:*n.  to  provide  f  lipfilr- 
in  areas  of  high  tauridar..  t-erm  incidence.  Thunders to;"?  are;,  however, 
difficult  to  identify  from  aDcv«=  in  tile  daytime  and  hence,  the  aircraft  ,v,} 
have  been  ebove  thunderstorms  ..n  occasions  xhtn  toe  pilots  were  unable  to 
identify  the  underlying,  clouds  as  tbur.derr  tonus.  Pieo,  tte  HICAT  program 
lacked  a  copsistert,  format  for  piloid  repvti.tq  thundarst.oiwr  when  in  smooth 
flight-.  / 

Clear  air  turbulence  veo  observed  in  the  atra'.  .'sphere  at  times  when  the: 
wer  ‘  no  tbunderot  errma  below  the  aircraft.  This  suggests  that  clear  air 
turbulence  observed  above  thunderstorms  is  not  necessarily  d.irectly  physic - 
clly  related  to  the  existence  of  a  thunderstorm  below.  The  record  of  'll 
observations  of  turbulence  on  of  42  reports  02"  thunders  toras  below  uhe 
aircraft  suggests  that  the  turbulence  and  the  thunderstorm  activity  J.3  signi¬ 
ficantly  related. 

World-Wide  Distribution  cf  Thu-ideratogms 

Thunderstorms  are  relatively  easily  identified  in  a  qualitative  etnae  vt  a 
quantitative  definition  is  difficult  device.  Dennis  (T;  defined  a 
thunderstorm  as  a  shower  complex  in  which  there  a one  or  .tore  cells  pro¬ 
ducing  lightning  and  within  wnieh  the  edge-to-edge  .:ell  separation  doe”  not 
exceed  3°  ka  (19  miles).  Some  area?  of  the  world  have  t  'Ofcdex arenas  con¬ 
sisting  of  only  one  cell  vLij.e  other  parts  of  the  v:  rlo  have  •  to.  t 
squall  lines  with  ««ny  ee3_Ia  ..  d  e.„G*wid Ing  1 -•  -•'■•'-.v  ui  Kilometers 

in  length.  _ _  a  given  cloud  mass  may  be  in  different  stages  of 

development*  One  cell  may  be  near  the  end  of  its  life  cycle  at  the  time  a 
nev  cell  is  starting  to  grow  rapidly.  The  duration  of  thunderstorms  Is  also 
variable.  Prophet  (8)  has  -shown  that  the  duration  increases  with  increasing 
size  of  the  radar  echo  area  of  the  thunderstorm.  For  our  present  purposes. 
Brooks' (9)  estimate  of  the  average  duration  as  one  hour  will  be  accepted. 

The  world -wide  distribution  of  thunderstorms  has  been  given  by  Brooks  (9), 
World  Meteorological  Organization  (10),  and  the  Handbook  of  Geophysics  (ll) 
in  terms  of  "thunderstorm  days".  A  thunderstorm  day  is  defined  as  a  local 
calendar  day  on  which  thunder  Is  heard.  Schonland  (12)  estimated  that  one 
out  of  seven  thunderstorm  days  have  more  than  one  thunderstorm  within 
audible  distance.  In  this  report  It  will  be  assumed  that  one  out  of  seven 
thunderstorm  days  has  two  thunderstorms. 

I 

The  literature  relating  to  the  thunderstorms  contains  a  number  of  estimates 
of  the  average  size  of  thunderstorms  and  of  the  area  seen  by  one  observer. 

The  estimates  of  these  two  areas  are  usually  so  nearly  equal  that  In  this 
report  they  will  be  assumed  to  be  equal.  To  determine  the  relative  pro¬ 
portion  of  the  total  area  of  a  region  that  1b  covered  by  thunderstorms  it 
is  not  necessary  to  know  the  numerical  value  of  the  area  of  the  thunder- 
stonns  because,  if  n  is  the  number  of  thunderstorm  days  for  a  given  period 
of  timey  then,  for  any  mean  area  of  the  thunderstorm,  the  ratio,  F,  of  the 
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j.  '-H  •••■{'  .  i  {«‘  _ '  '  j-  ,  txoe  in  turbulence  given  in  Section  5. 

xurvner,  these  tall  thunderstorms  may  be  seen  for  great  distances  and  can 
usually  be  easily  avoided.  Hence,  for  the  present,  the  evidence  indicates 
that  clear  air  turbulence  above  thunderstorms  is  real  but  thunderstorm 
penetration  into  the  stratosphere  is  infrequent  enough  to  Justify  not 
categorizing  thunderstorm  turbulence  in  a  separate  class. 
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Section  Vn 


A 


WORLD-WIDE  DI3TRI0UTIO*  CF  TOPXRAPKIC  CATEGORIES 

In  Sectious  II,  III,  and  V  data  were  presented  that  support  the  hypothesis 
that  the  underlying  topography  haa  a  significant  effect  upon  the  probability 
density  cf  the  FM5  gust  velocity,  the  scale  length, and  the  proportion  of  tiae 
In  turbulence  for  the  altitude  interval  undsr  consideration  (45,000  to 
65,000  ft.).  The  categories  cf  topography  were  defined  In  terms  of  differ¬ 
ences  in  elevation  and  not  of  elevation.  Figures  19  to  26  show  the  world¬ 
wide  distribution  of  these  topographic  categories.  Rote  that  there  Is  a 
sign*,  fi ,*nnt  difference  between  these  maps  and  naps  that  Indicate  altitude 
only.  Plateaus  are  shown  as  flatl«>nd  even  though  they  nay  be  at  a 
relatively  great  elevation. 

The  effects  of  mountain  ridges  nay  extend  to  as  .Mich  as  50  niles  leeward 
depending  upon  tbe  wind  speed  and  direction  aa?  upon  the  stability  of  the 
atmosphere. 

The  distribution  of  the  earth's  surface  area  in  terms  of  topographic  cats- 
gcrie*  is  « 3  follovn : 

VM  *  l*  71* 

Flavland  23 jt 

Low  Mountains  3 $ 

High  Mountains  jf> 

In  the  HICAT  nights  approximately  26%  of  the  night  miles  were  flown  above 
mountains,  a  figure  substantially  greater  than  the  6jt  of  the  earth's  surface 
area  in  this  category.  It  appears  that  the  proportion  of  time  in  turbulence 
expected  from  several  nights  randomly  distributed  throughout  the  world 
would  be  less  than  the  figures  recommended  In  Section  V,  Table  XVI,  for  the 
altitude  and  season  categories. 
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Figure  20.  South  America  -  Relief  Differences 
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Figure  21.  Europe  -  Relief  Differences 
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Figure  24,  Africa  -  Relief  Differences 


25.  Southeast  Asia  -  Relief  Differences 


Figure  26a.  Japan  -  Relief  Differences 
Figure  26b.  Hawaii  -  Relief  Differences 
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SECTION  vm 


FREQUENCY  OF  OCCURRENCE  CF  'TURBULENCE  AS  A  FUNCTION  OF  HEADWIND  AND  TAILWIND 

Turbulence  has  been  shown  to  have  occurred  more  frequently  during  periods 
when  the  aircraft  heading  was  upwind  (Aahburn,  Waco,  and  Mitchell,  (3)).  A 
further  Investigation,  using  only  the  264  runs  with  wind  speed  >  20  kts, 
has  shown  a  westerly  wind  component  In  86$  of  the  runs  (Figure  27A)  and  an 
aircraft  heading  with  a  westerly  component  6l$  of  the  time  In  turbulence 
(Figure  27B).  These  findings  would  Indicate  a  higher  expectancy  of  turbu¬ 
lence  with  a  heading  Into  the  wind.  This  is  shown  to  be  the  case  In  Figure 
27C  where  62$  of  the  flight  time  In  turbulence  coincided  with  a  headwind 
component.  This  Increased  to  67$  for  turbulence  >  moderate.  The  classifi¬ 
cation  of  turbulence  was  based  on  peaks  In  the  eg  acceleration  traces  (see 
Section  V). 

During  several  flights,  the  aircraft  was  flown  In  a  pattern  course,  en¬ 
countering  turbulence  more  than  once  In  the  same  area.  There  were  184  runs 
In  pattern  flights  where  the  wind  equalled  or  exceeded  20  kts.  Over  60$  of 
the  time  In  turbulence  was  again  associated  with  headwind  conditions.  This 
Indicates  that  time  of  day  or  selection  of  flight  path  cannot  be  used  to 
explain  the  Increased  frequency  of  occurrence  of  turbulence  with  a  headwind 
component. 

A  distinct  difference,  however,  occurs  between  the  percentage  of  time  In 
turbulence  during  headwind  conditions  for  the  early  flights  54  to  179  (67$) 
and  the  redirected  flights  180  to  285  (53 $).  Because  the  uncertainties  In 
deriving  winds  were  greater  in  the  early  phase  of  the  HIEAT  program,  it  is 
concluded  that  the  headwind-turbulence  relationship  may  be  explained,  in 
part,  as  a  result  of  errors  In  the  cr jputed  winds  in  the  first  Bet  of 
flights. 
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Figure  2JA.  Frequency  of  Occurrence  of  264-  Runs  by  Wind  Direction;  27B  Percent 
of  Time  in  Turbulence  by  Aircraft  Heading;  and  2JC  Percent  of  Time 
in  Turbulence  by  Aircraft  Heading  -  Wind  Direction  Difference. 
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SECTION  IX  • 


POWER  SPECTRAL  DENSITY  CURVES  AS  FUNCTIONS 
OF  AIRCRAFT  HEADING  WITH  RESPECT  TO  WIND  DIRECTION 

Introduction 

Several  spectra  with  long  wavelength  cutoffs  of  10,000  ft  (10"^  cycles/ft), 
normalized  by  passing  them  through  the  same  point  at  the  shortest  wavelength, 
were  compared  for  six  flights  (Figure  28).  Of  special  interest  was  the 
heading  orientation  with  respect  to  wind  direction.  Each  flight  contained 
spectra  from  runs  obtained  in  pattern  flight  over  the  same  general  area. 

Sets  of  two  flights  are  presented  for  high  mountains,  low  mountains,  and 
flatland  topographic  categories. 

High  Mountains 

(1)  Flight  114:  Figure  28A  shows  three  spectra  from  runs  east  of  the 
Sierra  Nevada  ridge  line.  Turbulence  was  light.  Winds  recorded  at  Oakland, 
California,  and  Wlnnemucca,  Nevada,  were  80  kts  from  the  southwest  at  20,000 
ft.  Runs  10  and  11  were  both  crosswind  and  nearly  overlapping  in  space,  yet 
produced  spectra  with  significantly  varying  slopes  for  \  >  2,000  ft.  All 
three  spectra  bend  upwards  with  slopes  from  1.67  to  2.64  at  longer  wave¬ 
lengths  compared  to  slopes  all  close  to  1.50  for  X  <  2,000  ft. 

(2)  Flight  280 :  Figure  28B.  The  three  runs  had  moderate  to  severe  turbulence, 
headinga  nearly  parallel  to  the  windflow  and  a  bending  upwards  (steeper  slopes) 
at  X  >  1,000  ft.  Run  10,  the  spectrum  with  the  largest  RMS  (2,000)  value,  1b 
the  least  steep  of  the  three.  The  flight  was  between  the  crest  of  the  Rockies 
and  Denver. 

Low  Mountains 

Flight  90 :  Figure  28C.  These  spectra  were  computed  from  data  obtained  over 
North  Island  in  New  Zealand.  Turbulence  was  light  to  moderate.  As  with  the 
previous  spectra  there  is  good  agreement  for  X  <  1,000  ft  but  variable  slopes 
at  long  wavelengths.  Again,  the  spectrum  with  the  highest  RMS  (2,000),  Run 
12, has  the  shallowest  slope  at  long  wavelengths.  Thfe  steepest  spectrum  (Run 
5)  has  the  lowest  RMS  (2,000). 

Flight  102:  Figure  28D.  These  two  spectra  agree  over  all  wavelengths  and 
represent  conditions  on  the  lee  side  of  the  Great  Dividing  Range  in  southeast 
Australia.  Moderate  turbulence  was  present.  Headings  appear  crosswind  in 
both  cases,  although  aircraft  winds  were  too  variable  to  establish  a  definite 
heading-wind  orientation.  The  spectra  bend  considerably  downward  beyond  X  * 
1,000  ft,  having  slopes  of  1.36  compared  to  1.64  at  short  wavelengths.  A 
similarity  exists  between  these  and  the  shape  of  the  spectra  for  Flight  90, 

Run  12.  The  runs  of  Flight  102  and  90,  Run  12  were  above  low  mountains 
(elevation  differences  <  4,000  ft)  whereas  the  other  samples  were  associated 
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with  mountains  having  relief  differences  >  6,000  ft. 

Flatland 

Flight  107:  Figure  28E.  Four  spectra  are  reproduced  here  from  data 
obtained  over  southeast  Australia  during  light  to  moderate  turbulent 
conditions.  Slopes  vary  appreciably  for  X  >  1,000  ft.  The  spectrum  from 
Run  6,  with  a  crosswind  component  in  the  heading,  behaves  similar  to  Run 
8  which  had  a  heading  nearly  parallel  to  the  wind.  In  contrast  to  the 
mountain  samples  Run  9,  with  the  highest  intensity,  has  the  steepest 
slope.  The  spectra  from  Runs  6  and  8  bend  considerably  for  X  >  1,000  ft. 
with  slopes  of  1.05  compared  to  1.66  for  X  <  1,000  ft. 

Flight  198:  Figure  28F.  These  were  from  samplings  over  thunderstorms  in 
the  southeast  United  States.  Slopes  vary  from  1.19  to  1.44  for  X  >  1,000 
ft.  and  are  close  to  I.60  at  short  wavelengths  for  all  spectra.  The  bending 
at  long  wavelengths  is  similar  in  magnitude  to  many  of  the  flatland  and  low 
mountain  samples.  All  of  Flight  198' s  spectra  were  crosswind  and  fairly 
high  in  intensity  showing  that  the  variation  in  shape  under  similar  orient¬ 
ations  can  be  large. 


Spectral  shapes  for  mountain  wave  samples  appear  to  depend  more  on  turbulence 
intensity  than  flight  orientation  with  respect  to  the  wind.  Low  intensity 
spectra  tend  to  bend  upward  at  long  wavelengths.  The  low  signal  to  noise 
ratio  at  long  wavelengths  could  be  a  factor  in  causing  the  upward  bending. 

The  variation  in  spectral  slopes  is  as  great  for  similar  headings  &b  for 
different  headings  with  respect  to  wind  direction.  The  shape  of  flatland 
spectra  also  show  no  correlation  between  aircraft  heading  and  wind  direction. 
Their  slopes  vary  considerably  at  long  wavelengths  and  in  general  are  much 


Figure  28A  and  28b.  Power  Spectral  Density  Curves  for  Various  Aircraft 

Headings  -  High  Mountains. 
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Figure  28C  and  28D.  Power  Spectral  Density  Curves  for  Various  Aircraft 

Headings  -  Low  Mountains. 
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SECTION  X 


DISTRIBUTION  OF  TRUE  GUST  VELOCITIES 


Introduction 


On  15  February  1 968  the  personnel  of  the  HICAT  Project  and  the  National 
Center  for  Atmospheric  Research  cooperated  in  an  investigation  of  turbu¬ 
lence  associated  with  mountain  waves  over  the  eastern  Rocky  Mountain  area 
near  Boulder,  Colorado.  Four  aircraft  participated  in  the  experiment. 

Each  aircraft  flew  a  track  between  Kremling  and  Akron,  Colorado.  Each 
aircraft  flew  at  different  altitudes.  The  HICAT  aircraft  flew  at  altitudes 
of  54,000  to  66,900  ft.  For  this  particular  flight  the  true  gust  velocities 
at  0.08  second  intervals  have  been  tabulated.  Trends  were  not  removed. 

Distribution  of  the  True  Gust  Velocities  for  HICAT  Flight  280,  Run  11. 

Flight  280,  Run  11  wgs  made  at  61,500  ft.  altitude.  The  heading  was  273° 
and  the  wind  was  276°  at  17  knots  (almost  a  direct  head  wind).  The  maxi¬ 
mum  U.  values  were  +  10. 5  and  -10. 3  ft/sec.  Root -mean -square  gust  velo- 
citiesefor  the  spectra  truncated  at  X  =  40,000  ft  (2.5  x  10 cycles/ft) 
are  4.21,  5*22,  and  5.48  ft/sec.  for  the  vertical,  lateral  and  longitudinal 
components,  respectively.  The  corresponding  root -mean-square  gust  velocities 
obtained  from  the  time  histories  with  trends  removed  were  4.87,  12.19  and 
10.86  ft/sec.  The  true  airspeed  was  419  knots  and  the  length  of  the  turbu¬ 
lent  period  was  539  seconds.  A  subjective  classification  of  light  to 
moderate  was  assigned  to  the  turbulence. 

Figure  29  illustrates  the  cumulative  frequencies  of  gust  velocities  of  given 
magnitudes  for  each  of  the  three  components.  Trends  were  not  removed.  If 
the  gust  velocities  were  normally  distributed  the  curves  would  be  straight 
lines,  ^he  curves  indicate  that  the  gust  velocities  are  not  normally 
distributed  and  that  the  turbulence  was  not  isotropic.  Computed  means  are 
1.19,  3»65,  and  -  0.00  ft/sec.  for  the  vertical,  lateral  and  longitudinal 
components,  respectively,  and  the  corresponding  standard  deviations  4.53;  12.6 
and  15.1  ft/sec.  The  data  used  to  construct  Figure  29  are  from  one  sample  of 
turbulence.  Analysis  of  other  data  is  required  before  broad  conclusions  may 
be  reached  on  the  distribution  of  gust  velocities. 

Crooks  et  al  (2)  list  the  root -mean -square  gust  velocities  computed  from 
the  time  histories  as  4.87,  12.19  and  10.86  ft/sec  for  the  vertical,  lateral 
and  longitudinal  components,  respectively.  These  are  significantly  different 
than  the  values  given  above.  The  difference  is  associated  with  data  sampling 
techniques,  removal  of  trends,  and  the  grouping  of  the  data  into  histograms. 
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Figure  29*  Distribution  of  True  Gust  Velocities  for  Three  Components  - 
Flight  280,  Run  11. 
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SECTION  XI 


SUlMUff  AND  CONCLUSIONS 

The  principal  goal  of  the  analyses  presented  in  this  report  was  to  determine 
numerical  values  for  the  standard  deviation  of  the  root -mean-square  gust 
velocities  (b.. ),  the  scale  lengths  and  slopes  of  the  pover  spectral  density 
curves,  and  trie  ratios  of  turbulent  flight  miles  to  total  flight  miles. 

These  numerical  values  are  given  in  Tables  XIH  and  XVI.  The  values  of 
are  shown  to  increase  significantly  with  roughness  of  terrain  for  all  three 
components.  Numerical  values  are  0.9  to  3*4  ft/sec  for  the  vertical  com¬ 
ponent,  1.5  to  4.6  ft/sec  for  the  lateral  component  and  1.2  to  4.6  ft/sec 

for  the  longitudinal,  component. 

There  was  no  well  defined  trend  in  the  variation  of  b.  with  altitude.  For 

the  total  H1CAT  sample  b.  was  equal  to  2.3,  3*0  and  3^0  ft/sec  for  the  verti¬ 

cal,  lateral  and  longitudinal  components,  respectively.  The  data  were  in¬ 
adequate  to  determine  if  a  further  breakdown  of  each  category  into  "storm" 
and  "nan -storm"  sets  would  be  useful.  The  values  of  the  b's  given  in  this 
report  are  dependent,  to  a  large  degree,  upon  the  choice  or  slopes,  m,  for 
the  power  spectral  density  curves.  When  the  b's  given  in  Table  XIII  are 
used,  it  is  important  to  use  the  appropriate  values  of  x  along  with  the  values 
of  b^. 

Composite  spectra  were  obtained  by  the  use  of  several  techniques,  each  of 
which  yield  different  values  of  m.  None  of  the  techniques  indicated  con¬ 
sistent  slopes  of  -5/3.  The  slopes  that  are  recommended  in  this  report 
(Table  XIII)  vary  from  -I.25  for  flights  over  water  to  -1.55  for  flights 
over  high  mountains  for  the  vertical  component  of  the  gust  velocity.  Com¬ 
parable  values  for  the  lateral  component  are  -I.50  and  -1-75  and  for  the 
longitudinal  component  are  -1.45  and. -1.60.  No  significant  change  of  m  with 
altitude  was  noted.  The  values  of  m  for  the  total  BICAT  sample  are  -1.45, 

-I.65  and  -1-55  for  the  vertical,  lateral  and  longitudinal  components,  res¬ 
pectively.  All  composite  spectra  indicated  relatively  large  scale  lengths. 

A  scale  length  of  at  least  4,000  ft.  for  all  three  components  Is  recommended 
as  best  fitting  the  data. 

Other  conclusions  presented  In  this  report  may  be  briefly  sumaarlzed  as  follows: 

1.  The  computed  values  of  the  P*  but  not  b*e  are  significantly  affected  by 
the  selection  of  the  lower  bound  of  turbulence. 

2.  The  ratio  of  turbulent  to  total  flight  miles  was  found  to  be  0.027  for 
flights  over  water  and  flatland,  0.036  for  flights  over  low  mountains 
and  0.049  for  flights  over  high  mountains  when  effects  due  to  pattern 
flying  were  removed.  These  ratios  were  reduced  approximately  35#  after 
-vami  ni  ng  figures  from  NASA  U-2  programs  which  contained  less  bias 
towards  turbulence  searching  techniques.  The  overall  ratio  of  turbu¬ 
lent  to  total  flight  miles  was,  in  addition,  considered  to  be  high 
because  of  the  larger  amount  of  miles  flown  over  mountains  (26#)  than 
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would  be  expected  on  a  mod  oh,  world -wide  basis  (6£). 

r 

3.  TUrbuleoce  decreased  substantially  with  altitude  over  flat  terrain 
but  showed  little  dropoff  with  altitude  over  Mountains.  Vinter 
appeared  to  be  the  season  with  naxlMua  turbulence. 

4.  Clear  air  turbulence  existed  above  thunderstorms  at  Infrequent  periods. 
Therefore,  no  separate  classification  was  Made  of  turbulence  assoc¬ 
iated  with  thunders tors  activity. 

5.  Root -Mean-square  gust  velocities  obtained  by  truncating  the  spectra 
at  X  *  2,000  ft  nay  In  Most  cases  be  assumed  to  have  a  normal  prob¬ 
ability  density  distribution,  except  at  large  values  of  the  EM3. 

6.  The  probability  density  distribution  of  the  Individual  gust  velocity 
values  (saapled  every  0.08  seconds)  has  been  shown  to  differ  from  a 
normal  distribution  In  one  case  over  Mountains. 

7.  Power  spectral  density  curves  for  selected  cases  over  various  topo¬ 
graphic  categories  showed  a  vide  range  of  shapes.  There  was  no 
apparent  correlation  between  spectral  shape  and  orientation  of  the 
aircraft  vlth  the  wind.  The  Mountain  wave  associated  spectra  had 
the  steepest  slopes  vlth  the  lower  Intensity  cases  showing  greater 
bending  upward. 


APFEKDIX  I 


DETERKDIATIOW  OF  SCALE  LEMCTH  BY  H»  RATIOS 


The  equation: 


#  (X)  =  (1  ♦  2  tt/X)  5/3 


when  integrated  between  a  truncated  X  and  the  shortest  X  ,  yields  the  RMS 
for  the  truncated  X  :  l/2 

>b 
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For  the  HBCAT  data,  X^  -  130  ft. 


An  expression  for  the  scale  length  obtained  by  solving  for  L  In  above  can 
be  rathei\  involved .  Consequently,  an  alternative  method  may  be  used  to 
arrive  atfL  which  Involves  calculating  ratios  of  selected  RMS  values  from 
""the  data  and  equating  these  to  similar  ratios  obtained  from  the  family  of 
mild  knee  equations. 
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where  X  was  selected  as  4,000  ft.  or  10,000  ft.  and  ^  j  qqq  ^ 

*i  2  "■  i/s 

1  •  1  _ 

(  1  +  0.000628LJ2/3  (  1  +  0.04ii3L}2/3 
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X  =  10,000  ft) 
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T  i'  +  0*00&28l773  (  1 "+  0.0^31)2/3 
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IS.  ABSTRACT 

Report  consists  of  analysis  of  measurements  made  in  HICAT  flight  program  and  describee 
in  reports  Issued  in  1967  and  1968.  Clear  air  turbulent  regions  in  the  altitude  range 
45,000  to  65,000  ft.  varied  in  length  from  1  to  115  nm.  Power  spectral  density  curvet 
computed  from  data  from  the  longer  turbulent  regions  had  slopes  and  scale  lengths  thal 
varied  over  a  relatively  vide  range.  Composite  spectra  for  altitude  and  topographic 
categories  for  the  long  turbulent  regions  had  slopes  of  -1.25  to  -1.55  for  the  verti¬ 
cal  component,  -I.50  to  -I.70  for  the  lateral  component,  and  -1.45  to  -I.60  for  the 
longitudinal  component.  The  scale  lengths  were  all  above  8,000  ft.  The  standard 
deviations  of  the  rcot -mean-square  gust  velocities  increase  slgnlglcantly  with  in¬ 
crease  in  roughness  of  underlying  terrain.  The  range  in  values  is  0.9  to  3.4  ft/sec. 
for  the  vertical  component,  1.5  to  4.6  ft/sec  for  the  lateral  component  and  1.2  to 
4.5  ft/sec.  for  the  longitudinal  component.  There  was  relatively  little  change  with 
altitude.  The  minimum  root -mean-square  velocity  to  be  considered  turbulence  varied 
O.95  to  3*02  ft/sec.  if  truncated  spectra  were  used  to  establish  basic  definitions  of 
turbulence. 
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